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ABSTRACT 

We present the results of a determination of the galaxy luminosity func- 
tion at ultraviolet wavelengths at redshifts of z — 0.0 — 0.1 from GALEX 
data. We determined the luminosity function in the GALEX FUV and NUV 
bands from a sample of galaxies with UV magnitudes between 17 and 20 that 
are drawn from a total of 56.73 deg 2 of GALEX fields overlapping the bj- 
selected 2dF Galaxy Redshift Survey. The resulting luminosity functions are 
fainter than previous UV estimates and result in total UV luminosity densi- 
ties of 10 25 - 55±0 - 12 ergs s- 1 Hz" 1 Mpc~ 3 and 10 25 - 72±0 - 12 ergs s" 1 Hz" 1 Mpc" 3 at 
1530A and 2310A, respectively. This corresponds to a local star formation rate 
density in agreement with previous estimates made with Ha-selected data for 
reasonable assumptions about the UV extinction. 

Subject headings: surveys - galaxies: luminosity function - ultraviolet: galaxies 
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1. Introduction 

In the past few years determinations of the star formation history of the universe have 
allowed us to begin to understand quantitatively when and how the stars in the universe 
were formed. Measurements of the rest-frame ultraviolet luminosities of galaxies have 
been particularly useful in this endeavor. In the very local universe, there is a relative 
lack of systematic surveys of galaxies in the UV. Before the launch of GALEX, the most 
comprehensive survey of galaxies in the local universe was from the FOCA experiment 
(Milliard et al. 1992), a balloon-borne telescope that made measurements in a single band 
centered at 2000A. Based upon FOCA observations of a total of ~ 2.2 deg 2 , Treyer et 
al. (1998) and Sullivan et al. (2000) measured the first UV luminosity function (LF) for 
a sample of 273 galaxies with spectroscopic redshifts at z = 0.15. Their LF has a steep 
faint end slope and a total UV luminosity density, and corresponding star formation rate 
density, larger than most previous estimates. This higher local UV luminosity density in 
conjunction with measurements at larger distances lead Wilson et al. (2002) to infer a 
luminosity density evolution proportional to (1 + z) L7±1 "°, a trend shallower than had been 
estimated previously from the CFRS sample (Lilly et al. 1996). 

In this letter we present the first results regarding the UV LF based upon measurements 
from the Galaxy Evolution Explorer (GALEX) in conjunction with redshifts from the 2dF 
Galaxy Redshift Survey (2dFGRS) (Colless et al. 2001). The new GALEX data allow us to 
expand upon the previous FOCA results using a much larger sample drawn from an area 
of 56.73 deg 2 although to a shallower limiting magnitude of m uv = 20. Throughout this 
paper, we assume H = 70 km s -1 Mpc" 1 , Vl M = 0.3 and = 0.7. 
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2. Data 

The data analyzed in this paper consist of 133 GALEX All-Sky Survey (AIS) 
pointings that overlap the 2dF Galaxy Redshift Survey in the South Galactic Pole region. 
The GALEX field-of-view is circular with diameter of 1.2° and each pointing is imaged 
simultaneously in both the FUV and NUV bands with effective wavelengths of 1530A and 
2310A, respectively. The median exposure time for the fields is 105 seconds, allowing us to 
reach a S/N ratio of ~ 5 for FUV « 20.0 and NUV « 20.5. See Martin et al. (2004) and 
Morissey et al. (2004) for details regarding the GALEX instruments and mission. 

Sources were detected and measured from the GALEX images using the program 
SExtractor (Bertin & Arnouts 1996). As the NUV images are substantially deeper than 
the FUV, we used the NUV images for detection and measured the FUV flux in the same 
aperture as for the NUV. The fields analyzed here were processed using a larger SExtractor 
deblending parameter DEBLEND_MINCONT as the standard GALEX pipeline processing 
tends to break apart well-resolved galaxies into more than one source. We elected to use 
the MAG_AUTO magnitudes measured by SExtractor through an elliptical aperture whose 
semi-major axis is scaled to 2.5 times the first moment of the object's radial profile, as first 
suggested by Kron (1980). All of the apparent magnitudes were corrected for foreground 
extinction using the Schlegel et al. (1998) reddening maps and assuming the extinction law 
of Cardelli et al. (1989). The ratio of the extinction in the GALEX bands to the reddening 
E(B — V) was calculated by averaging the extinction law over each GALEX bandpass, 
resulting in A FUV /E(B - V) = 8.376 and A NUV /E(B - V) = 8.741. The median extinction 
correction for the galaxies in our South Galactic Pole sample is 0.15 mag in both bands, 
with the corrections ranging from 0.1 to 0.3 mag. 

The GALEX catalogs were matched with the 2dFGRS input catalog using a search 
radius of 6". To remove any overlap between adjacent pointings, we only included sources 
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detected within the inner 0.45° of each field. In addition, sources likely contaminated by 
artifacts from bright stars, with 2dF redshift quality flag less than three or with effective 
exposure times less than 60 sec were removed. Finally, we excluded GALEX sources in 
regions where the 2dF redshift completeness was less than 80%. After applying all of these 
cuts to each band, the total area on the sky of GALEX-2dF overlap is 56.73 deg 2 . 

The GALEX resolution of 6 - 7" (FWHM) (Morissey et al. 2004) is not sufficient to 
accurately separate stars and galaxies. Furthermore, the 2dFGRS input catalog available 
from the 2dFGRS web page 1 only includes galaxies brighter than bj = 19.45 and does not 
include stars. In order to asses the total completeness of our 2dF-GALEX matched sample, 
we normalized our results to the total galaxy number counts determined by Xu et al. (2004) 
based primarily upon 22.64 deg 2 of GALEX Medium Imaging Survey data overlapping the 
Sloan Digital Sky Survey (SDSS) Data Release 1 (Abazajian et al. 2003). As the SDSS 
data include stars and galaxies and reach fainter magnitudes, they result in a more accurate 
determination of the galaxy number counts in the UV. If we assume that the average galaxy 
number counts in the SDSS North Galactic Pole fields are the same as in the GALEX-2dF 
overlap, then the redshift completeness of the 2dF matched catalog is given by the number 
counts of galaxies with redshifts from 2dF divided by the total galaxy number counts from 
the SDSS overlap. This ratio is shown in Figure 1. 

The completeness turns over at 20th mag because the redshift sample becomes 
incomplete for galaxies with blue (FUV — bj) or (NUV — bj) colors. We have limited our 
LF determination in each band to galaxies brighter than this limit. To avoid problems with 
photometry of large bright galaxies, we also imposed a bright magnitude limit of 17. The 
average completeness weighted by the number counts in the range 17 — 20 mag is 92% in 
the FUV and 79% in the NUV. For objects with magnitudes brighter than 20.0 in either 

1 http://www. mso.anu.edu/2dFGRS/ 
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band, we visually inspected all of the 2dF spectra and removed a total of 27 objects with 
very broad emission lines indicating that the objects are most likely some sort of AGN. The 
redshift distributions for the FUV and NUV samples are shown in Figure 2. 

We further restricted our sample to those galaxies with redshifts z < 0.1 to insure that 
our sample is not sensitive to evolution. The average redshifts are 0.055 and 0.058 in the 
FUV and NUV, respectively. After applying all of the cuts mentioned in this section, a 
total of 896 galaxies in the FUV and 1124 galaxies in the NUV remained. The luminosity 
functions for the objects with z > 0.1 are presented in Treyer et al. (2004). 



Using the FUV, NUV and bj magnitudes, we assigned a best-fit spectral type to each 
galaxy using a representative subset of the SEDs from Bruzual & Chariot (2003) and 
determined the K-correction needed to transform the observed UV magnitudes to rest-frame 
measurements at z — 0. The K-corrections are in general quite small (< 0.2). 

We determined the LF $(M) and its error <t($(M)) in each band using the V max 
method (Felten 1976): 



where f(m) is the inverse of the redshift completeness as estimated in §2 above and V max is 
the maximum co-moving volume within which each galaxy could have been observed given 
the bright and faint limiting magnitudes of our sample and its best-fit SED. The resulting 
LFs are shown in Figure 3. 



By minimizing x 2 , we fit the V max LF points in each band with a Schechter function 
(Schechter 1976): §(L)dL = <p*(L/ L*) a e~ L / L * dL/ L* where <p\ M* and a were free 



3. Luminosity Functions 
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Fig. 1. — Completeness of the GALEX-2dF redshift sample denned as the ratio of the 
number counts of galaxies with 2dF redshifts to the number counts of galaxies as derived 
from GALEX observations that overlap the SDSS survey (Xu et al. 2004). The solid and 
dashed lines indicate the FUV and NUV redshift completeness, respectively. 
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Fig. 2. — The redshift distributions of the FUV and NUV selected samples (blue solid and 
red dashed lines, respectively) in the range 17 < m uv < 20. 
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parameters. The best fit parameters and their errors, calculated using the range of solutions 
within 1.0 of the minimum x 2 , are listed in Table 1 along with the best-fit LF from Sullivan 
et al. (2000) converted to the AB magnitude system and to H = 70. The errors in a and 
M* are highly correlated and the inset of Figure 3 shows the la error contours projected 
into the M* — a plane. Since the V max method can be biased in the presence of clustering, 
we also computed the best-fit Schechter parameters using the maximum likelihood STY 
method (Sandage et al. 1979). The resulting STY values are listed in Table 1 and are 
also plotted in the inset of Figure 3. The STY values lie just inside and outside of the la 
Vmax error ellipses in the FUV and NUV, respectively. We adopt the V max results in the 
discussion below. 



4. Discussion 

As can be seen in Figure 3, there are significant differences between the results 
presented here and those from Sullivan et al. (2000). The GALEX results have a fainter M* 
in both bands and have a shallower faint end slope. The FOCA passband is centered at 
2015A with FWHM of 188A and thus one would expect the FOCA results to lie in between 
the GALEX FUV and NUV data. However, the FOCA sample is truly an UV-selected 
sample while that presented here relies upon the ^-selected 2dFGRS. This selection could 
introduce a bias in our results if the galaxies for which we do not have redshifts have a 
different redshift distribution than the galaxies which are included in our sample. On the 
other hand, it is now well established that the UV luminosity density increases with redshift 
(e.g. Somerville et al. 2001) and part of the difference is likely a real effect (Treyer et al. 
2004). However, the difference of ~ 0.9 mag between the FOCA and NUV values for M* 
would require evolution much larger than determined from other surveys as well as GALEX 
data at higher redshifts (Arnouts et al. 2004; Schiminovich et al. 2004). A preliminary 
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Fig. 3. — The FUV (blue circles) and NUV (red triangles) LFs for z < 0.1. The solid lines 
are the Schechter function fits with best-fit parameters from Table 1. The dotted green line 
shows the LF measured at 2000 A from FOCA data by Sullivan et al. (2000) over the range 
of absolute magnitudes explored in that study. The inset plots the la error contours of the 
Schechter function fits projected into the M* — a plane for the FUV (blue) and NUV (red). 
The dashed contour shows values with x 2 ~ Xmin = 1-0 while the solid contour delineates 
X 2 — Xmin = 2-3 which corresponds to the joint la uncertainty on M* and a. The red and 
blue stars indicate the best-fitting values of M* and a obtained from the STY method. 



Table 1. Schechter Function Parameters 

V max method STY method 

band z M* a log0* (Mpc -3 ) logp L (ergs s" 1 Hz" 1 Mpc -3 ) M* a 

FUV 0-0.1 -18.04 ±0.11 -1.22 ±0.07 -2.37 ±0.06 25.55 ±0.12 -18.12 -1.23 

NUV 0-0.1 -18.23 ±0.11 -1.16 ±0.07 -2.26 ± 0.06 25.72 ±0.12 -18.27 -1.10 

FOCA 0.15 -19.10 ±0.13 -1.51 ±0.10 -2.48 ±0.11 26.06 ±0.15 
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comparison of the GALEX and FOCA photometry in a couple of overlapping fields indicates 
that the FOCA magnitudes are on average brighter with the difference becoming larger for 
fainter sources. It appears likely these offsets and non-linearities in the FOCA photometry 
account for a major part of the difference between the FOCA and GALEX LFs with the 
remainder likely due to a combination of galaxy evolution and the FOCA sample selection. 

In Table 1 we also list the total luminosity density calculated from the best-fit Schechter 
parameters as pi = J °° L$(L)dL = <p*L*Y{a + 2). The statistical errors in logpL that take 
into account the covariance between the three Schechter function parameters are 0.02 in 
each band. In addition to this error, the uncertainty in the GALEX photometric zeropoint 
is ft 10% in both bands, corresponding to an uncertainty in logp L of 0.04. A potentially 
larger source of error is that due to large scale structure. The variation in the number 
density of galaxies n in a contiguous volume V is given approximately by 5n/n ft (J3/V) 1 / 2 
(Davis & Huchra 1982) where J 3 is an integral over the galaxy 2-point correlation function 
and has a value of ~ 10 4 Mpc 3 for a correlation function of the form £(r) = (r/r )~ 7 with 
ro = 7.21Mpc and 7 = 1.67 (Hawkins et al. 2003). The galaxy number counts from Xu 
et al. (2004) used to set the normalization of our LFs were derived from approximately 
22.64 deg 2 . For z < 0.1, the corresponding rms variation in the number density would be 
5n/n ft 0.24, or an uncertainty in 5\ogp L ft 0.11. Since UV-selected, star-forming galaxies 
are likely less clustered than optically-selected samples, this value is really an upper limit. 
Adding these uncertainties due to large scale structure and calibration in quadrature to the 
statistical errors results in a total uncertainty of 5\ogpi ft 0.12 in both bands. 

The spectral slope (3, defined as fx oc X 13 with f\ in units of ergs s _1 A -1 Mpc -3 , 
corresponding to our two luminosity density measurements is (5 ft —1.1. This is slightly 
bluer than the slope of /3 — —0.9 determined by Cowie et al. (1999) at z — 0.7 — 1.3 from 
measurements at longer rest frame wavelengths spanning 1700A to 2750A. 
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The FUV luminosity density can be used to estimate the star formation rate (SFR) 
density in the local universe. For a constant star formation history and a Salpeter 
IMF, the SFR is related to the UV luminosity L v (in the range 1500 - 2800A) by 
SFR (M yr- 1 ) = 1.4 x 10 _28 L„(ergs s" 1 Hz" 1 ) (Kennicutt 1998). For the FUV luminosity 
density in Table 1, we obtain log (SFR FUV )(M yr" 1 Mpc~ 3 ) = -2.30 ± 0.12 with no 
extinction correction. For comparison, the extinction-corrected Ha LF at z < 0.045 
from Gallego et al. (1995) shifted to our assumed Hubble constant corresponds to 
log (SFRh«) = -1.86 ± 0.04 using the Ha to SFR conversion from Kennicutt (1998). 
Based upon Ha imaging of a subsample of the galaxies used by Gallego et al. (1995), 
Perez-Gonzalez et al. (2003) argued that the local Ha luminosity density is ~ 60% higher 
due to uncertainties in the aperture corrections applied to the spectroscopic data and 
corresponds to log (SFR Ha ) = — 1.6 ± 0.2. Bringing the FUV SFR into agreement with this 
result would require an extinction of A FUV — 1.8. 

An average extinction of A FUV ~ 1.8 is consistent with a simple estimate made using 
the observed (FUV — NUV) colors. While there is a well-defined relationship between the 
UV extinction and the spectral slope for starburst galaxies, more quiescent galaxies tend to 
have less extinction for a given UV slope than would be inferred from nearby starbursts 
(Bell 2002). In particular Kong et al. (2004) used the population synthesis models of 
Bruzual & Chariot (2003) along with the prescription described in Chariot & Fall (2000) for 
determining how starlight is absorbed by dust in galaxies to show that the smaller extinction 
in non-starburst galaxies can be explained by variations in the galaxies' star formation 
histories. Based upon a set of Monte Carlo realizations of these models spanning a range of 
extinctions, ages and star formation histories, Kong et al. (2004) were able to approximate 
the dependence of the FUV extinction on the UV spectral slope (3 with the following 
formula: A FUV = 3.87 + 1.87(/3 + 0.40 log b) where the variable b parametrizes the star 
formation history and is defined as the ratio of current to the past average star formation 
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rate. Assuming a constant star formation history (6 = 1), an extinction of Apuv = 1-8 is 
obtained for a spectral slope (5 = —1.1, a value consistent with that measured from the 
FUV and NUV luminosity densities. On the other hand, the average (FUV — NUV) color 
of our FUV-selected sample is 0.14, corresponding to a spectral slope of f3 — —1.67. For this 
P and 6=1, the Kong et al. (2004) formula results in A FUV = 0.7 mag. This extinction 
is similar to the results of Buat et al. (2004) who found that the average extinction for a 
local NUV-selected sample is A FUV ~ 1 mag based upon the FIR to UV flux ratio. If an 
extinction of A FUV ~ 1 is more appropriate for the UV-selected sample presented here, 
then the UV-based star formation rate density would be log (SFR FU v) = —1-9 ± 0.1, a 
value lower than that from Ha although still consistent to within the errors. In reality the 
extinction is likely a function of absolute magnitude and future GALEX papers will address 
in more detail correcting UV fluxes for extinction in a more rigorous way. 

In the near future we will continue our investigation of the UV luminosity function in 
the local universe using GALEX AIS data covering ~ 1000 deg 2 of the SDSS. In addition 
to expanding our sample to include more galaxies, we will use the SDSS photometry and 
spectroscopy to explore the dependence of UV luminosity on other galaxy characteristics, 
such as color, surface brightness, environment, metallicity and stellar mass. 

GALEX (Galaxy Evolution Explorer) is a NASA Small Explorer, launched in April 
2003. We gratefully acknowledge NASA's support for construction, operation, and science 
analysis for the GALEX mission, developed in cooperation with the Centre National 
d'Etudes Spatiales of France and the Korean Ministry of Science and Technology. 
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